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(57) Abstract: A method, system, and computer program product related to the diagnosis of diabetes, and is directed to predicting 
the long-term risk of hyperglycemia, and the long-term and short-term risks of severe hypoglycemia in diabetics, based on blood 
glucose readings collected by a self-monitoring blood glucose device. The method, system, and computer program product pertain 
directly to the enhancement of existing home blood glucose monitoring devices, by introducing an intelligent data interpretation 
component capable of predicting both HbAic and periods of increased risk of hypoglycemia, and to the enhancement of emerging 
continuous monitoring devices by the same features. With these predictions the diabetic can take steps to prevent the adverse conse- 
quences associated with hyperglycemia and hypoglycemia. 
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METHOD, SYSTEM, AND COMPUTER PROGRAM PRODUCT 
FOR THE EVALUATION OF GLYCEMIC CONTROL IN 
DIABETES FROM SELF-MONITORING DATA 

5 

CROSS-REFERENCES TO RELATED APPLICATIONS 

The present invention claims priority from U.S. Provisional Patent Application 
Serial No. 60/193,037 filed March 29, 2000, entitled "Algorithm for the Evaluation of 
10 Glycemic Control in Diabetes From Self-Monitoring Data" the entire disclosure of 
which is hereby incorporated by reference herein. 

US GOVERNMENT RIGHTS 

This invention was made with United States Government support under Grant 
15 Nos. NIH / NIDDK: ROl DK 28288 and NIH / NIDDK: ROl DK 51562, both 

awarded by National Institutes of Health. The United States Government has certain 
rights in the invention. 

FIELD OF THE INVENTION 
20 The present system relates generally to Glycemic Control of individuals with 

diabetes, and more particularly to a computer-based system and method for evaluation 
of predicting glycosylated hemoglobin (HbAic and HbAi) and risk of incurring 
hypoglycemia. 

25 BACKGROUND OF THE INVENTION 

Extensive studies, including the Diabetes Control and Complications Trial 
(DCCT) {See DCCT Research Group: The Effect Of Intensive Treatment Of Diabetes 
On The Development And Progression Of Long-Term Complications Of Insulin- 
Dependent Diabetes Mellitus. New England Journal of Medicine, 329: 978-986, 

30 1993), the Stockholm Diabetes Intervention Study (See Reichard P, Phil M: Mortality 
and Treatment Side Effects During Long-term Intensified Conventional Insulin 
Treatment in the Stockholm Diabetes Intervention Study. Diabetes, 43: 313-317, 
1994), and the United Kingdom Prospective Diabetes Study {See UK Prospective 
Diabetes Study Group: Effect of Intensive Blood Glucose Control With Metformin On 
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Complications In Patients With Type 2 Diabetes (UKPDS 34). Lancet, 352: 837-853, 
1998), have repeatedly demonstrated that the most effective way to prevent the long 
term complications of diabetes is by strictly maintaining blood glucose (BG) levels 
within a normal range using intensive insulin therapy. 
5 However, the same studies have also documented some adverse effects of 

intensive insulin therapy, the most acute of which is the increased risk of frequent 
severe hypoglycemia (SH), a condition defined as an episode of neuroglycopenia 
which precludes self-treatment and requires extemal help for recovery {See DCCT 
Research Group: Epidemiology of Severe Hypoglycemia In The Diabetes Control and 

10 Complications Trial. American Journal of Medicine, 90: 450-459, 1 991 , and DCCT 
Research Group: Hypoglycemia in the Diabetes Control and Complications Trial. 
Diabetes, 46: 271-286, 1997). Since SH can result in accidents, coma, and even death, 
patients and health care providers are discouraged from pxirsuing intensive therapy. 
Consequently, hypoglycemia has been identified as a major barrier to improved 

15 glycemic control (Cryer PE: Hypoglycemia is the Limiting Factor in the Management 
Of Diabetes. Diabetes Metab Res Rev, 15: 42-46, 1999). 

Thus, patients wdth diabetes face a life-long optimization problem of 
maintaining strict glycemic control without increasing their risk of hypoglycemia. A 
major challenge related to this problem is the creation of simple and reliable methods 

20 that are capable of evaluating both patients' glycemic control and their risk of 
hypoglycemia, and that can be applied in their everyday environments. 

It has been well known for more than twenty years that glycosylated 
hemoglobin is a marker for the glycemic control of individuals with Diabetes Mellitus 
(Type I or Type II). Numerous researchers have investigated this relationship and 

25 have found that glycosylated hemoglobin generally reflects the average BG levels of a 
patient over the previous two months. Since in the majority of patients with diabetes 
the BG levels fluctuate considerably over time, it was suggested that the real 
connection between integrated glucose control and HbAic would be observed only in 
patients known to be in stable glucose control over a long period of time. 

30 Early studies of such patients produced an almost deterministic relationship 

between the average BG level in the preceding 5 weeks and HbAic, and this curvilinear 
association yielded a correlation coefficient of 0.98 (See Aaby Svendsen P, Lauritzen 
T, Soegard U, Nerup J (1982). Glycosylated Hemoglobin and Steady-State Mean 
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Blood Glucose Concentration in Type 1 (Insulin-Dependent) Diabetes, Diabetologia , 
23, 403-405). In 1993 the DCCT concluded that HbAic was the "logical nominee" for 
a gold-standard glycosylated hemoglobin assay, and the DCCT established a linear 
relationship between the preceding mean BG and HbAic {See Santiago JV (1993). 
5 Lessons from the Diabetes Control and Complications Trial, Diabetes , 42, 1549-1554). 
Guidelines were developed indicating that an HbAic of 7% corresponds to a 
mean BG of 8.3 mM (150 mg/dl), an HbAic of 9% corresponds to a mean BG of 11. 7 
mM (210 mg/dl), and a 1% increase in HbAic corresponds to an increase in mean BG 
of 1 .7 mM (30 mg/dl, 2). The DCCT also suggested that because measuring the mean 
10 BG directly is not practical, one could assess a patient's glycemic control with a single, 
simple test, namely HbAic. However, studies clearly demonstrate that HbAic is not 
sensitive to hypoglycemia. 

Indeed, there is no reliable predictor of a patient's immediate risk of SH from 
any data. The DCCT concluded that only about 8% of future SH could be predicted 
15 from known variables such as the history of SH, low HbAic, and hypoglycemia 

unawareness. One recent review details the current clinical status of this problem, and 
provides options for preventing SH, that are available to patients and their health care 
providers (See BoUi, GB: How To Ameliorate The Problem of Hypoglycemia In 
Intensive As Well As Nonintensive Treatment Of Type I Diabetes. Diabetes Care, 22, 
20 Supplement 2: B43-B52, 1999). 

Contemporary home BG monitors provide the means for frequent BG 
measurements through Self-Monitoring of BG (SMBG). However, the problem with 
SMBG is that there is a missing link between the data collected by the BG monitors, 
and HbAic and hypoglycemia. In other words, there are currently no reliable methods 
25 for evaluating HbAic and recognizing imminent hypoglycemia based on SMBG 
readings {See Bremer T and Gough DA: Is blood glucose predictable from previous 
values? A solicitation for data. Diabetes 48:445-451, 1999). 

Thus, an object of this invention is to provide this missing link by proposing 
three distinct, but compatible, algorithms for evaluating HbAic and the risk of 
30 hypoglycemia from SMBG data, to be used to predict the short-term and long-term 
risks of hypoglycemia, and the long-term risk of hyperglycemia. 

The inventors have previously reported that one reason for a missing link 
between the routinely available SMBG data and the evaluation of HbAic and the risk 



3 



wo 01/72208 PCT/USOl/09884 
of hypoglycemia, is that the sophisticated methods of data collection and clinical 
assessment used in diabetes research, are infrequently supported by diabetes-specific 
and mathematically sophisticated statistical procedures. 

Responding to the need for statistical analyses that take into account the 
5 specific distribution of BG data, the inventors developed a symmetrizing 

transformation of the blood glucose measurement scale {See Kovatchev BP, Cox DJ, 
Gonder-Frederick LA and WL Clarke (1997). Symmetization of the Blood Glucose 
Measurement Scale and Its Applications, Diabetes Care , 20, 1655-1658) that works as 
the follows. The BG levels are measured in mg/dl in the United States, and in mmol/L 
10 (or mM) in most other coxintries. The two scales are directly related by 18 mg/dl = 1 
mM. The entire BG range is given in most references as 1.1 to 33.3 mM, and this is 
considered to cover practically all observed values. According to the 
recommendations of the DCCT (See DCCT Research Group (1993) The Effect Of 
Intensive Treatment of Diabetes On the Development and Progression of Long-Term 
1 5 Complications of Insulin-Dependent Diabetes Mellitus. New England Journal of 
Medicine^ 329, pp 978-986) the target BG range - also known as the euglycemic 
range- for a person v^th diabetes is 3.9 to 10 mM, hypoglycemia occurs when the BG 
falls below 3.9 mM, and hyperglycemia is when the BG rises above 10 mM. 
Unfortunately, this scale is numerically asymmetric — the hyperglycemic range (10 to 
20 33.3mM) is wider than the hypoglycemic range (LI to 3.9mM), and the euglycemic 
range (3.9 to lOmM) is not centered within the scale. The inventors correct this 
asymmetry by introducing a transformation, /(^5G>, which is a continuous function 
defined on the BG range [L 1, 33. 3J, having the two-parameter analytical form: 

/(BG, a, P) = [(In (BG)f -pl a, P> 0 
25 and which satisfies the assumptions: 

Al : f(333. a, p) = -/(7.i, a, p) and 
A2: /r/aft a, P) = -/(3.P, a, P), 
Next,/(C^ is multiplied by a third scaling parameter to fix the minimum and 

maximum values of the transformed BG range at - VlO and VlO respectively. These 
30 values are convenient since a random variable with a standard normal distribution has 
99.8% of its values within the interval [- VlO , VlO ] . If BG is measured in mmol/1, 
when solved numerically with respect to the assumptions Al and A2, the parameters of 
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the function f(BG, a, JS) eacQ a= 1, 026, fi = 1.861, and the scaling parameter is r = 
1 . 794, If BG is measured in mg/dl instead, the parameters are computed to be « = 
L084,yff-5.557, and 1,509, 

Thus, when BG is measured in mmoi/l, the symmetrizing transformation is 
5 f(BG) = L 794[(ln f^G;/^^^- L861]. and when BG is measured in mg/dl the 
symmetrizing transformation isf(BG) = 1.509 [(In (BG))^'^^"* - 5,381], 

On the basis of the symmetrizing transformation f(,) the inventors introduced 
the Low BG Index - a new measure for assessing the risk of h3^oglycemia from 
SMBG readings (See Cox DJ, Kovatchev BP, Julian DM, Gonder-Frederick LA, 
1 0 Polonsky WH, Schlundt DG, Clarke WL: Frequency of Severe Hypoglycemia In 
IDDM Can Be Predicted From Self-Monitoring Blood Glucose Data. Journal of 
Clinical Endocrinology and Metabolism, 79: 1659-1662, 1994, and Kovatchev BP, 
Cox DJ, Gonder-Frederick LA Young-Hyman D, Schlundt D, Clarke WL. Assessment 
of Risk for Severe Hypoglycemia Among Adults With IDDM: Validation of the Low 
15 Blood Glucose Index, Diabetes Care 21:1870-1875, 1998). Given a series of SMBG 
data the Low BG Index is computed as the average of 10, f(BGf taken for values of 
f(BG) <0 and 0 otherwise. Also suggested was a High BG Index, computed in a 
symmetrical to the Low BG Index manner, however this index did not find its practical 
application. 

20 Using the Low BG Index in a regression model the inventors were able to 

accoimt for 40% of the variance of SH episodes in the subsequent 6 months based on 
the SH history and SMBG data, and later to enhance this prediction to 46% (See 
Kovatchev BP, Straume M, Farhi LS, Cox DJ: Estimating the Speed of Blood Glucose 
Transitions and its Relationship With Severe Hypoglycemia. Diabetes, 48: 

25 Supplement 1, A363, 1999). 

In addition, the inventors developed some data regarding HbAic and SMBG 
(See Kovatchev BP, Cox DJ, Straume M, Farhy LS. Association of Self-monitoring 
Blood Glucose Profiles with Glycosylated Hemoglobin. In: Methods in Enzvmology. 
vol. 321 : Numerical Computer Methods, Part C, Michael Johnson and Ludvig Brand, 

30 Eds., Academic Press, NY; 2000). 

These developments became a part of the theoretical background of this 
invention. In order to bring this theory into practice, several key theoretical 
components, among other things, as described in the foUowdng sections, were added. 
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In particular, three methods were developed for employing the evaluation of HbAic, 
long-term and short-term risk for hypoglycemia. The development of these methods 
was, but not limited thereto, based on detailed analysis of data for 867 individuals with 
diabetes that included more than 300,000 SMBG readings, records of severe 
5 hypoglycemia and determinations of HbAi c. 

The inventors have therefore sought to improve upon the aforementioned 
limitations associated with the conventional methods, and thereby provide simple and 
reliable methods that are capable of evaluating both patients' glycemic control and 
their risk of hypoglycemia, and that can be applied in their everyday environments. 

10 

SUMMARY OF THE INVENTION 

The invention includes a data analysis method and computer-based system for 
the simultaneous evalviation, from routinely collected SMBG data, of the two most 
important components of glycemic control in diabetes: HbAic and the risk of 

15 hypoglycemia. For the purposes of this docxmient, self-monitoring of BG (SMBG) is 
defined as any method for determination of blood glucose at diabetic patients' natural 
environment and includes the methods used by contemporary SMBG devices 
customarily storing 200-250 BG readings, as well as methods used by emerging 
continuous monitoring technologies. Given this broad definition of SMBG, this 

20 invention pertains directly to the enhancement of existing home blood glucose 

monitoring devices by introducing an intelligent data interpretation component capable 
of predicting both HbAic and periods of increased risk of hypoglycemia, as well as to 
enhancement of future continuous monitoring devices by the same features. 

One aspect of the invention includes a method, system, and computer program 

25 product for evaluating HbAu from a predetermined period of collected SMBG data, 
for example 4-6 weeks. In one embodiment, the invention provides a computerized 
method and system for evaluating the HbAic of a patient based on BG data collected 
over a predetermined duration. The method includes computing weighted deviation 
toward high blood glucose (WR) and estimated rate of change of blood glucose (Dr) 

30 based on the collected BG data; estimating HbAic using a predetermined mathematical 
formula based on the computed WR and Dr; and providing a predetermined confidence 
interval for classification of said estimated value of HbAic. 
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Another aspect of the invention includes a method, system, and computer 
program product for estimating the long-term probability for severe hypoglycemia 
(SH). This method uses SMBG readings from a predetermined period, for example 4- 
6 weeks, and predicts the risk of SH within the following 6 months. In one 
5 embodiment, the invention provides a computerized method and system for evaluating 
the long term probability for severe hypoglycemia (SH) of a patient based on BG data 
collected over a predetermined duration. The method includes: computing weighted 
deviation toward low blood glucose (WL) and estimated rate of fall of blood glucose in 
the low BG range (DrDn) based on the collected BG data; estimating the number of 

10 future SH episodes using a predetermined mathematical formula based on the 

computed WL and DrDn; and defining a probability of incurring a select number of 
SH episodes respective to said estimated SH episodes. 

Still yet another aspect of the invention includes a method, system, and 
computer program product for identifying 24-hour periods (or other select periods) of 

15 increased risk of hypoglycemia. This is accomplished through the computation of the 
short-term risk of hypoglycemia using SMBG readings collected over the previous 24 
hours. In one embodiment, the invention provides a computerized method and system 
for evaluating the short term risk for severe hypoglycemia (SH) of a patient based on 
BG data collected over a predetermined duration. The method includes: computing 

20 weighted deviation toward low blood glucose (WL); determining Max(wl) by 
calculating maximum value of wl(BG;2)\ determining risk value by taking the 
geometric mean of WL and Max(wl) over the predetermined duration; providing a 
predetermined threshold risk value; and comparing the determined risk value to the 
threshold risk value. 

25 These three aspects of the invention can be integrated together to provide 

continuous information about the glycemic control of an individual with diabetes, and 
enhanced monitoring of the risk of hypoglycemia. 

These and other objects, along with advantages and features of the invention 
disclosed herein, will be made more apparent from the description, drawings and 

30 claims that follow. 
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BRIEF DESCRIPTION OF THE DRAWINGS 

The foregoing and other objects, features and advantages of the present 
invention, as well as the invention itself, will be more fully understood from the 
following description of preferred embodiments, when read together with the 
5 accompanying drawings in which: 

FIG. 1 is a flow chart illustrating the method of calculating the estimated HbAic 
and predicted HbAic confidence intervals in accordance with the present invention. 

FIG. 2 is a flow chart illustrating the method of calculating the estimated 
number of fiiture SH episodes and the associated probability thereof in accordance 
1 0 with the present invention. 

FIG. 3 is a flow chart illustrating the method of calculating the estimated short 
term risk of an incurring imminent SH in accordance with the present invention. 

FIG. 4 is graphical representation of a typical BG disturbances observed before 
and after an episode of severe hypoglycemia. 
15 FIG. 5 illustrates the action of the method for predicting short-term SH by 

presenting 10 weeks of data for Subject A (upper panel) and Subject B (lower panel). 
SH episodes are marked by triangle; a black line presents the risk value. When the risk 
threshold is crossed, the method indicates a subsequent high-risk period (gray bar). 

FIG. 6 is a functional block diagram for a computer system for implementation 
20 of the present invention. 

FIGS. 7-9 are schematic block diagrams of altemative variations of the present 
invention related processors, communication links, and systems. 

DETAILED DESCRIPTION OF THE INVENTION 

25 The invention makes possible, but not limited thereto, the creation of precise 

methods for the evaluation of diabetics' glycemic control, and include, firmware and 
software code to be used in computing the key components of the method. The 
inventive methods for evaluating HbAu, the long-term probability of SH, and the 
short-term risk of hypoglycemia, are also validated based on the extensive data 

30 collected, as will be discussed later in this document. Finally, the aspects of these 
methods can be combined in structured display or matrix. 
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Stationary measures of BG deviation 
According to the inventors' theory of BG symmetrization (See Kovatchev BP, 
Straume M, Cox DJ, Farhi LS. Risk Analysis of Blood Glucose Data: A Quantitative 
Approach to Optimizing the Control of Insulin Dependent Diabetes. J of Theoretical 
5 Medicine , 3:1-10,^2001) the natural clinical center of the BG measurement scale is at 
BG level of 1 12.5 mg/dl (6.25 mmol/1) - a safe euglycemic value for a diabetes 
patient. 

Given this clinical center of the BG scale, the weighted deviations to the left 
(towards hypoglycemia) or to the right (towards hyperglycemia) are computed. The 
10 degree of weighting of these deviations will be represented by parameters a and b 
respectively as follows: 

wl(BG;a) = lO.ffBG/ if f(BG)<0 and 0 otherwise, and 
wr(BG:b)^ 10.f(BGf if /fBG;>0 and 0 otherwise, 
where/(^5G^ is the BG symmetrization function presented in the background section. 
15 The weighting parameter a and b could be different, or the same for the left and right 
deviations. The inventors* data analyses demonstrated that the optimal for practical 
application parameter values are a=2 (which is the parameter value used for 
computation of the Low BG Index) and 6=7. Given a series of BG readings xu xz 
Xn, the average weighted deviations to the left and to the right of the clinical center of 
20 the BG scale are defined as: 

WL = w/(x,;2) and WR = iv>vr(x,;l) respectively. 

w /=i 

These two measures of BG deviation do not depend on the timing of the BG readings, 
and therefore are stationary. In order to capture the dynamics of BG change, measures 
of the BG rate of change are introduced as provided below. 

25 

Computation of BG risk rate of change 
Letxy, ... x„ be « SMBG readings of a subject recorded at time points 
tu ... tn^ This data is next transformed by calculating the numbers f(xi), /(x^) , 
f( Xn) and draw a cubic spline S(t) passing through the points (ti J(xi)), (t2j(x2,)) 
30 ( tn J( x„)). Thus, the fimction S(t) is a continuous fimction defined on the whole 
interval [t] , trj and such that S(tj)=f(xj), for Also calculated are the set of 
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numbers Sk =10,S(k+ tjf for MJ, r„ - thus getting interpolated values at one- 
hour increments. 

Next, consider all couples of numbers s^ with consecutive indices: 

Co=(so,sj), C]=(sj,s^ , C2^(szS3X and denote by M^^the set of all couples 

5 such that Sk > Sk-\-i and by the set of all couples such that Sk < Sk^i . 

Finally, let DrDn be the average of the numbers Sk-^i - Sk, provided that Q e 
Mdn , and Dr be the average of the numbers s^+i - Sk, provided that Q e Mup+ Mdn^ 

The numbers DrDn and Dr provide a measure for the rate of change of BG in a 
"risk space," e.g. the rate of change of the risk associated with any BG level change. 
10 In addition, DrDn measures the rate of BG change only when BG goes down, i.e. 
DrDn evaluates how quickly the risk could increase when BG falls, while Dr is a 
measure of the overall vulnerability of BG to fluctuations. It is further asserted that 
DrDn will be associated with risk for hypoglycemia (if someone's blood glucose could 
fall quickly, his/her risk for hypoglycemia would be higher), while Dr will be 
15 associated with the overall stability of BG. 

Software code Tpresented in SPSS control language^ 
The first is for when the BG readings are in mmol/L, and in this case the 
variable is BGMM. The second is for when the BG readings are in mg/dl, and in this 
20 case the variable is BGMG. 

If BG is measured in mmol/L, each BG reading is first transformed as follows: 

SCALEl=(ln(BGMM) )**1.026 - 1.861 
RISK1=32 . 185*SCALE1*SCALE1 

If BG is measured in mg/dl, each BG reading is first transformed as follows: 

SCALE2= (In (BGMG) ) **1. 08405 - 5.381 
RISK2=22 . 765*SCALE2*SCALE2 

Further, the left and right weighted deviations are computed as follows: 

WL=OWL=0 

IF (SCALEl le 0.0) WL=RISK1 
30 WR=0 

IF (SCALEl gt 0.0) WR=sqrt (RISKl ) 
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Provided that the BG readings are equally spaced in time, or are interpolated at 
one-hour increments, the BG rate of change is computed as: 

Dr = RISK1{BG)-RISK1{BG-1) 
DrDn=0 

5 IF (SCALE le 0.0 and Dr gt 0) DrDn=Dr 

Finally, an aggregation pass through all BG readings for a subject will produce: 

WL = mean(WL) 
WR = mean{WH) 

Dr = mean(Dr), and DrDn = mean(DrDn) 

10 

Method for the Evaluation of HbA i^ 
A preferred embodiment of HbAic evaluation method 100 according to the 
invention is illustrated in FIG. 1 . In a first step 1 02, SMBG data is collected over a 
predetermined period of time. For example, the SMBG data is collected over 4-6 

15 weeks with a frequency of 3-5 BG measurements per day, of which are transformed by 
the code or formulas presented in the previous section. Different formulas are to be 
used if the BG measurements are stored in mg/dl, or in mmol/l. One skilled in the art 
would appreciate that various levels, durations, and frequencies can be employed. In a 
step 104, weighted deviation towards high blood glucose (WR) and estimated rate of 

20 change of blood glucose (Dr) is computed using the formula / code discussed above. 
In a step 106, an estimate of HbAic from self-monitoring data is computed using the 
linear fimction: EstHBAlc = 0.9008* WR - 0.8207*DR + 6.7489. It is noted that the 
coefficients of this fimction are derived from data for 867 individuals with diabetes, 
and one would recognize that fiirther data accumulation may update these coefficients. 

25 In step 108 HbAic estimate categories representing a range of values for estimated 
HbAic are defined according to Table 1. 



Table 1: Defining categories on the basis of EstHBAlc: 



EstHBAlc 


<7.8 


7.8 - 8.5 


8.5 - 9.0 


9.0 - 9.6 


9.6-10.3 


10.3-11.0 


> 11.0 


Category 


1 


2 


3 


4 


5 


6 


7 
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In Step 110 predicted confidence intervals for corresponding HbAic are derived 
according to Table 2. 



Table 2: Predicted 95% confidence intervals for classification of HbAici 



Category 


1 


2 


3 


4 


5 


6 


7 


HBAic 


<8.0 


8.0-8.5 


8.5-9.0 


9.0 - 9.5 


9.5-10.1 


10.1-11.0 


> 11.0 



5 In step 1 12, the estimated HbAic fi^om step 106 is assigned in one of the 

categories provided in Table 1 and / or Table 2. 

Empirical Validation of Evaluation of HbA ir 
The intervals for HbAic in Table 2 are based on extensive research. To validate 

10 these intervals we analyzed SMBG and HbAic data from 867 subjects v^ith diabetes. 
All subjects were instructed to use BG memory meters for six months and to measure 
their BG two to four times a day. During the same period 5 to 8 HbAic assays were 
performed for each subject. The memory meter data were electronically downloaded 
and stored in a computer for fiirther analysis. This procedure produced a database 

15 containing more than 300,000 SMBG readings and 4,180 HbAic assays taken over six 
months. Analysis of variance was conducted to compare HbAic in the seven categories 
identified in Table 1 . The five categories were highly significantly different, with 
F=91 and p<0.00001. Moreover, the average HbAic was significantly different for 
each pair of categories as demonstrated by Duncan's ranges, with p<0.01 . 

20 Also, 95% confidence intervals were computed for the mean value of HbAic in 

each of the seven categories. These confidence intervals were used as a basis for 
computing the HbAic intervals presented in Table 2. Post-hoc analysis of the 
classification power of this method demonstrated that the method was well protected 
against extreme errors such as incorrectly classifying HbAic in category 1, 2 or 3 on 

25 the basis of SMBG while the actual HbAic was greater than 9.5%, or classifying HbAic 
in category 5, 6 or 7 while the actual HbAic was below 9.0%. 

In summary, after an initial 4-6 weeks of SMBG readings the computerized 
method computes an interval estimate for the value of HbAic that can be used to track 
patients* changes in glycemic control in the high BG range. 
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Method for Evaluation of the Long-Term Probability for Severe Hypoglycemia (SIT) 
A preferred embodiment of long-term probability for SH evaluation method 
200 according to the invention is illustrated in FIG. 2. In a first step 202, SMBG data 
is collected over a predetermined period of time. For example, the SMBG data is 
5 collected over 4-6 weeks with a frequency of 3-5 BG measurements per day, of which 
are transformed by the code or formulas presented immediately above. Different 
formulas are to be used if the BG measurements are stored in mg/dl, or in mmol/l; 
One skilled in the art would appreciate that various levels, durations, and frequencies 
can be employed. In a step 204, WL and DrDn are computed using the formula / code 
10 as discussed above. In step 206, an estimate of the nvimber of future SH episodes is 
computed using the linear function: 

EstNSH = 3.3613*WL - 4.3427*DrDn - 1.2716. 
It is noted that the coefficients of this function are derived fi"om data for 1 81 
individuals with diabetes, and one would appreciate that further data accimiulation 
15 may update these coefficients. It is further noted that this formula provides a single 
value estimate for the number of future SH episodes and that through additional 
methodologies, as discussed below, categories are provided with ranges and 
confidence levels for enhanced clinical applications. In step 208, estimated number of 
SH episodes (estNSH) categories representing a range of values for estNSH are 
20 defined according to Table 3. 

Table 3: Classification of EstNSH: 



EstNSH 


< 0.775 


0.775 - 3.750 


3.750-7.000 


> 7.000 


Category 


1 


2 


3 


4 



25 



In step 210, respective to the estNSH categories, the probability of incurring 0, 
1 -2, or more than 2 SH episodes in the following six months is derived, as represented 
in table 4. 





Category 1 


Category 2 


Category 3 


Category 4 


0 SH 


90% 


50% 


25% 


< 20% 


1-2 SH 


10% 


25% 


25% 


>2SH 


25% 


50% 


> 80% 
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In step 212, the EstNSH from step 206 is assigned in one of the categories 
provided in Table 3 and / or Table 4. 



Empirical Validation of Evaluation of the Long-Term Probability for SH 
5 One-hundred-eighty-one adults with Type 1 diabetes (mean age 37 years, 

duration of diabetes 18 years) used memory meters to collect more than 34,000 SMBG 
over a month. Then for the next six months they recorded in diaries any occxirrence of 
SH. The SMBG data were mathematically transformed and an a linear regression 
model was used to predict future severe hypoglycemia resulting in a highly significant 

10 model (F=36.3, p<0.0001) and multiple R of 55%, 

All subjects were classified into 4 categories using the present long-term SH 
method. The average number of future SH episodes in categories 1, 2, 3, and 4 was 
0.3, 2.0, 5.0, and 9.75 respectively. Analysis of variance demonstrated highly 
significant differences between these categories, F=19.0, p<0.0001. 

15 In summary, a linear combination of the Low BG Index and the rate of drop of 

BG as measured in "risk space" provide an accurate assessment of the long-term risk 
of SH. Because it is based on SMBG records that are automatically stored by many 
reflectance meters, this is an effective and clinically useful indicator of patients' 
glycemic control in the low BG range. 

20 

Method for the Evaluation of the Short-term (within 24 hours) Risk of Hypoglycemia 

A preferred embodiment of short term risk of SH evaluation method 300 
according to the invention is illustrated in FIG. 3. In a first step 302, SMBG is data is 
collected over a predetermined short term period. For example, the SMBG data is 

25 collected over a 24 hour period, with a frequency of 3-5 BG measurements per day - 4 
or more readings, as a nominal level according to data analyses. One skilled in the art 
would appreciate that various levels, periods (durations), and frequencies can be 
employed. In a step 304 WL(24) and Max(wl) is computed from all readings collected 
within the preceding 24 hours, wherein the maximum value of wl(BG;2) is Max(wl). 

30 In step 306, the risk value is by taking the geometric mean of WL and Max(wl) over 
the 24 hour period, wherein said risk value is mathematically defined as: 

Risk(24) = ^WL(24)^Max{wl)\ 



14 



wo 01/72208 PCT/USOl/09884 
In step 308 a threshold risk value is determined. In step 3 10 the estimated risk 
value is compared to the threshold risk value. For example, if the threshold risk value 
is set at 17, then if Risk(24)>l 7, then-based on the SMBG data collected over the 
previous 24 hours- the resultant indication is a high risk of the patient incurring 
5 imminent hypoglycemia. In other words, this is a decision-making rule that considers 
a 24-hour period of SMBG data and judges whether this period is likely to precede an 
imminent hypoglycemia episode. The threshold value of 17 is derived from an 
extensive data set, however, it is recognized that it is possible that this value maybe 
adjusted with further accumulation of data or for additional objectives. 

10 

Empirical Validation of Evaluation of the Short-term Risk of Hvpoglvcemia: 
Eighty-five individuals were recruited through advertisement in newsletters, 
diabetes clinics, and through direct referrals. The inclusion criteria were: 1) age of 21- 
60 years; 2) type I diabetes with at least two years duration, and insulin use since the 

15 time of diagnosis; 3) at least 2 documented SH episodes in the past year; and 4) routine 
use of SMBG devices for diabetes monitoring. The participants were instructed to use 
the meter 3-5 times a day, and to record in monthly diaries any SH episodes, including 
the exact dates and times of their occurrences. SH was defined as severe 
neuroglycopenia that results in stupor or imconsciousness and precludes self-treatment. 

20 For each subject the study continued 6-8 months and each month the subject's meter 
was downloaded and the SH diary was collected. The memory capacity of the meters 
was sufficient, and the downloading was often enough, so that no BG data were lost. 
No changes were made in the participants' diabetes management routine, nor were any 
additional treatments administered during the study. 

25 During the study a total of 75,495 SMBG readings (on average 4.0±1 .5 per 

subject per day) were downloaded from the participants' memory meters, and 399 
(4.7±6.0 per subject) SH episodes were recorded in their diaries. An important 
finding, among other things, was that episodes of moderate or severe hypoglycemia are 
preceded and followed by measurable BG disturbances. In the 24-hour period before 

30 an SH episode the Low BG Index (e.g. WL) rose (p<0.001), the average BG was lower 
(p=0.001), and the BG variance increased (p=0.001). In the 24 hours following the SH 
episode, the Low BG Index and BG variance remained elevated (p<0.001), but the 
average BG returned to its baseline. 
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To this end, FIG. 4 is graphical representation of a typical BG disturbance 
observed before and after an episode of severe hypoglycemia. In the period 48 to 24 
hours before the SH episode, the average BG level decreased and the variance of BG 
increased. In the 24-hour period immediately preceding the SH episode, the average 
5 BG level dropped further and the variance of BG continued to increase. In the 24-hour 
period following the SH episode, the average BG level normalized, but the BG 
variance remained greatly increased. Both the average BG and its variance retumed to 
their baseline levels v^ithin 48 hours after the SH episode. 

As such, as part of the invention, the disturbances presented in FIG. 4 are 

10 quantified from SMBG data to enable the evaluation of the short-term risk of 

hypoglycemia. The cutoff value of Risk(24)=l 7 is derived from an optimization along 
the following restrictions: 1) the method had to predict a maximum percentage of SH 
episodes, i.e. to identify as risky a maximum percentage of 24-hour periods preceding 
SH, and 2) to prevent overestimation of the risk, the method had to identify as risky no 

15 more that 15% of the total time of the study (one day a week on average). The cutoff 
risk value of 17 was held constant for all subjects. The reason for choosing the value 
of 1 5% was to prevent the patients from becoming irritated with an overabundance of 
"false alarms" and then ignoring ''true alarms." In practice, a patient's physician can 
select an altemate value depending on the severity of the patient's diabetes and 

20 particular objectives. 

The following example illustrates the action of the algorithm on the SMBG 
data of two participants in the study. FIG. 5 presents ten weeks of data for Subject A 
(upper panel) and Subject B (lower panel). SH episodes are marked by triangles; a 
black curve presents the risk value. When the risk threshold (the horizontal line at 

25 Risk=l 7) is crossed, the algorithm indicates a subsequent high-risk period (gray bar). 
For Subject A, 7 out of 9 SH episodes are predicted and there are 5 false alarms, e.g. 
high-risk periods that did not result in SH; for Subject B there are 3 false alarms and 
the only SH episode is predicted. It is obvious that Subject B's risk values when 
compared to Subject A's risk values, include more and higher deviations. For both 

30 subjects, all SH episodes were accompanied by supercritical risk values, and about half 
of all large deviations were accompanied by one or more SH episode. 

Across all participants in the study, 44% of all recorded SH episodes were 
preceded, within 24 hours, by a high-risk period, and 50% were preceded, within 48 
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hours, by a high-risk period. If only periods with either at least 3, or at least 4 SMBG 
measurements were considered, the accviracy of the latter prediction increased to 53% 
and 57%, respectively. Post-hoc analysis of BG levels occurring during, or 
immediately after, high-risk periods that were not followed by an SH episode, i.e. 
5 during, or immediately after false alarms, demonstrated that the average per subject 
minimum of such BG levels was 2.3±0.2 mmol/1 versus 5.9±L7 mmol/1 (t=19.5, 
p<0.0001) for all non-risk periods, including all SH episodes that remained 
unaccounted for. This indicates that, although symptomatic SH did not occur, BG 
levels following high-risk periods were notably low. 
10 In summary, the inventors simulated the action of the short-term risk method 

on a 6-month series of SMBG readings for 85 individuals with Type I diabetes. With 
four or more SMBG readings per day, at least 50% of all episodes of SH could be 
anticipated. Even when symptomatic SH did not occur, the algorithm predicted 
episodes of moderate hypoglycemia. 

15 

Integration of the Three Methods 
The three methods of this invention, as discussed above and illustrated in FIGS. 
1-3, utilize the same series of SMBG data. Therefore, from an SMBG-device point of 
view, a unified display or matrix of the results of these three methods could be made 
20 similar to the grid output presented below: 





EstHBA categories (Algorithm 1) 


1 


2 


3 


4 


5 


6 


7 


EstNSH 
categories 


1 








i: Ss 1-. . - 
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Thus, for example, the output for subject 1 (Ss 1) shown in the above grid 
indicates that this person is likely to have HbAic between 9 and 9.5%, and has a 90% 
25 chance not to experience severe hypoglycemia in the subsequent 6 months. The output 
for subject 2 (Ss 2) indicates that this person is likely to have HbAic below 8%, and 
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has a greater than 80% chance to experience at least 3 SH episodes in the subsequent 6 
months. 

In addition to this grid-output, the short term risk method provides a continuous 
tracking of the risk of imminent hypoglycemia and can be used to sound an alarm 
5 when this risk becomes high. 

The method of the invention may be implemented using hardware, software or 
a combination thereof and may be implemented in one or more computer systems or 
other processing systems, such as personal digit assistants (PDAs). In an example 
embodiment, the invention was implemented in software running on a general ptirpose 

10 computer 900 as illustrated in FIG. 6, Computer system 600 includes one or more 
processors, such as processor 604. Processor 604 is connected to a communication 
infi-astructure 606 (e.g., a commimications bus, cross-over bar, or network). Computer 
system 600 includes a display interface 602 that forwards graphics, text, and other data 
from the communication infrastructure 606 (or from a frame buffer not shown) for 

1 5 display on the display unit 630. 

Computer system 600 also includes a main memory 608, preferably random 
access memory (RAM), and may also include a secondary memory 610. The 
secondary memory 610 may include, for example, a hard disk drive 612 and/or a 
removable storage drive 614, representing a floppy disk drive, a magnetic tape drive, 

20 an optical disk drive, etc. The removable storage drive 614 reads from and/or writes to 
a removable storage unit 618 in a well known manner. Removable storage unit 618, 
represents a floppy disk, magnetic tape, optical disk, etc. which is read by and written 
to by removable storage drive 614. As will be appreciated, the removable storage unit 
618 includes a computer usable storage mediimi having stored therein computer 

25 software and/or data. 

In alternative embodiments, second£iry memory 610 may include other means 
for allowing computer programs or other instructions to be loaded into computer 
system 600. Such means may include, for example, a removable storage unit 622 and 
an interface 620. Examples of such removable storage xmits/interfaces include a 

30 program cartridge and cartridge interface (such as that foimd in video game devices), a 
removable memory chip (such as a ROM, PROM, EPROM or EEPROM) and 
associated socket, and other removable storage xmits 622 and interfaces 620 which 
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allow software and data to be transferred from the removable storage unit 622 to 
computer system 600. 

Computer system 600 may also include a communications interface 624. 
Communications interface 624 allows software and data to be transferred between 
5 computer system 600 and extemal devices. Examples of communications interface 624 
may include a modem, a network interface (such as an Ethemet card), a 
communications port, a PCMCIA slot and card, etc. Software and data transferred via 
communications interface 624 are in the form of signals 628 which may be electronic, 
electromagnetic, optical or other signals capable of being received by commimications 
10 interface 624. Signals 628 are provided to commimications interface 624 via a 

communications path (i.e., channel) 626. Channel 626 carries signals 628 and may be 
implemented using wire or cable, fiber optics, a phone line, a cellular phone link, an 
RF link and other communications channels. 

In this docimaent, the terms "computer program medium" and "computer usable 
15 medium" are used to generally refer to media such as removable storage drive 914, a 
hard disk installed in hard disk drive 612, and signals 628. These computer program 
products are means for providing software to computer system 600. The invention 
includes such computer program products. 

Computer programs (also called computer control logic) are stored in main 
20 memory 608 and/or secondary memory 610. Computer programs may also be received 
via communications interface 624. Such computer programs, when executed, enable 
computer system 600 to perform the features of the present invention as discussed 
herein. In particular, the computer programs, when executed, enable processor 604 to 
perform the functions of the present invention. Accordingly, such computer programs 
25 represent controllers of computer system 600. 

In an embodiment where the invention is implemented using software, the 
software may be stored in a computer program product and loaded into computer 
system 600 using removable storage drive 614, hard drive 612 or commimications 
interface 624. The control logic (software), when executed by the processor 604, 
30 causes the processor 604 to perform the fimctions of the invention as described herein. 

In another embodiment, the invention is implemented primarily in hardware 
using, for example, hardware components such as application specific integrated 
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circuits (ASICs). Implementation of the hardware state machine to perform the 
functions described herein will be apparent to persons skilled in the relevant art(s). 

In yet another embodiment, the invention is implemented using a combination 
of both hardware and software. 
5 In an example software embodiment of the invention, the methods described 

above were implemented in SPSS control language, but could be implemented in other 
programs such as, but not limited to, C + + programming language. 

FIGS. 7-9 show block diagrammatic representation of altemative 
embodiments of the invention. Referring FIG. 7, there is shown a block diagrammatic 

10 representation of the system 710 essentially comprises the glucose meter 728 used by a 
patient 712 for recording, inter alia, insulin dosage readings and measured blood 
glucose ("BG") levels. Data obtained by the glucose meter 728 is preferably 
transferred through appropriate communication links 714 or data modem 732 to a 
processing station or chip, such as a personal computer 740, PDA, or cellular 

15 telephone. For instance, data stored may be stored within the glucose meter 728 and 
may be directly downloaded into the personal computer 740 through an appropriate 
interface cable. An example is the ONE TOUCH monitoring system or meter by 
LifeScan, Inc. which is compatible with IN TOUCH software which includes an 
interface cable to down load the data to a personal computer. 

20 The glucose meter is common in the industry and includes essentially any 

device that can functions as a BG acquisition mechanism. The BG meter or 
acquisition mechanism, device, tool, or system includes various conventional methods 
directed toward drawing a blood sample (e.g. by fingerprick) for each test, and a 
determination of the glucose level using an instrument that reads glucose 

25 concentrations by electromechanical or claorimetric methods. Recently, various 
methods for determining the concentration of blood analytes without drawing blood 
have been developed. For example, U.S. Pat. No. 5,267,152 to Yang et al. describes a 
noninvasive technique of measuring blood glucose concentration using near-IR 
radiation diffuse-reflection laser spectroscopy. Similar near-IR spectrometric devices 

30 are also described in U.S. Pat. No. 5,086,229 to Rosenthal et al. and U.S. Pat. No. 
4,975,581 to Robinson et al. 

U.S. Pat. No. 5,139,023 to Stanley describes a transdermal blood glucose 
monitoring apparatus that relies on a permeability enhancer (e.g., a bile salt) to 
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facilitate transdermal movement of glucose along a concentration gradient established 
between interstitial fluid and a receiving medium. U.S. Pat. No. 5,036,861 to 
Sembrowich describes a passive glucose monitor that collects perspiration through a 
skin patch, where a cholinergic agent is used to stimulate perspiration secretion from 
5 the eccrine sweat gland. Similar perspiration collection devices are described in U.S. 
Pat. No. 5,076,273 to Schoendorfer and U.S. Pat. No. 5,140,985 to Schroeder. 

In addition, U.S. Pat. No. 5,279,543 to Glikfeld describes the use of 
iontophoresis to noninvasively sample a substance through skin into a receptacle on 
the skin surface. Glikfeld teaches that this sampling procedure can be coupled with a 
10 glucose-specific biosensor or glucose-specific electrodes in order to monitor blood 
glucose. Moreover, Intemational Publication No. WO 96/001 10 to Tamada describes 
an iontophoretic apparatus for transdermal monitoring of a target substance, wherein 
an iontophoretic electrode is used to move an analj^e into a collection reservoir and a 
biosensor is used to detect the target analyte present in the reservoir. Finally, U.S. Pat. 
1 5 No. 6,144,869 to Bemer describes a sampling system for measuring the concentration 
of an analyte present. 

Further yet, the BG meter or acquisition mechanism may include indwelling 
catheters and subcutaneous tissue fluid sampling. 

The computer or PDA 740 includes the software and hardware necessary to 
20 process, analyze and interpret the self-recorded diabetes patient data in accordance 
with predefined flow sequences (as described above in detail) and generate an 
appropriate data interpretation output. Preferably, the results of the data analysis and 
interpretation performed upon the stored patient data by the computer 740 are 
displayed in the form of a paper report generated through a printer associated with the 
25 personal computer 740. Altematively, the results of the data interpretation procedure 
may be directly displayed on a video display imit associated with the computer 740. 

FIG. 8 shows a block diagrammatic representation of an alternative 
embodiment having a diabetes management system that is a patient-operated apparatus 
810 having a housing preferably sufficiently compact to enable apparatus 810 to be 
30 hand-held and carried by a patient. A strip guide for receiving a blood glucose test strip 
(not shown) is located on a sxirface of housing 816. Test strip is for receiving a blood 
sample from the patient 812. The apparatus includes a microprocessor 822 and a 
memory 824 connected to microprocessor 822. Microprocessor 22 is designed to 
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execute a computer program stored in memory 824 to perform the various calculations 
and control functions as discussed in great detail above. A keypad 816 is connected to 
microprocessor 822 through a standard keypad decoder 826. Display 814 is connected 
to microprocessor 822 through a display driver 830. Microprocessor 822 
5 communicates vdth display driver 830 via an interface, and display driver 830 updates 
and refreshes display 814 under the control of microprocessor 822. Speaker 854 and a 
clock 856 are also connected to microprocessor 822. Speaker 854 operates under the 
control of microprocessor 822 to emit audible tones alerting the patient to possible 
future hypoglycemia. Clock 856 supplies the current date and time to microprocessor 
10 822. 

Memory 824 also stores blood glucose values of the patient 812, the insulin 
dose values, the insulin types, and the parameter values used by microprocessor 822 to 
calculate futvire blood glucose values, supplemental insulin doses, and carbohydrate 
supplements. Each blood glucose value and insulin dose value is stored in memory 

15 824 with a corresponding date and time. Memory 824 is preferably a non-volatile 
memory, such as an electrically erasable read only memory (EEPROM). 

Apparatus 810 also includes a blood glucose meter 828 connected to 
microprocessor 822. Glucose meter 828 is designed to measure blood samples received 
on blood glucose test strips and to produce blood glucose values from measurements of 

20 the blood samples. As mentioned previously, such glucose meters are well known in 
the art. Glucose meter 828 is preferably of the type which produces digital values 
which are output directly to microprocessor 822. Altematively, blood glucose meter 
828 may be of the type which produces analog values. In this altemative embodiment, 
blood glucose meter 828 is connected to microprocessor 822 through an analog to 

25 digital converter (not shown). 

Apparatus 810 further includes an input/output port 834, preferably a serial 
port, which is connected to microprocessor 822. Port 834 is connected to a modem 832 
by an interface, preferably a standard RS232 interface. Modem 832 is for establishing 
a communication link between apparatus 810 and a personal computer 840 or a 

30 healthcare provider computer 838 through a communication network 836. Specific 
techniques for connecting electronic devices through connection cords are well known 
in the art. Another altemative example is "bluetooth" technology communication. 
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Alternatively, FIG. 9 shows a block diagrammatic representation of an 
alternative embodiment having a diabetes management system that is a patient- 
operated apparatus 910, similar as shown in FIG. 8, having a housing preferably 
sufficiently compact to enable the apparatus 910 to be hand-held and carried by a 
5 patient. However, the present embodiment includes a separate or detachable glucose 
meter or BG acquisition mechanism 928. 

Accordingly, the embodiments described herein are capable of being 
implemented over data commimication networks such as the intemet, making 
evaluations, estimates, and information accessible to any processor or computer at any 
1 0 remote location, as depicted in FIGS, 6-9 and/or U.S. Pat. No. 5,85 1 , 1 86 to Wood, of 
which is hereby incorporated by reference herein. Alternatively, patients located at 
remote locations may have the BG data transmitted to a central healthcare provider or 
residence, or a different remote location. 

In summary, the invention proposes a data analysis computerized method and 
15 system for the simultaneous evaluation of the two most important components of 
glycemic control in individuals with diabetes: HbAic and the risk of hypoglycemia. 
The method, while using only routine SMBG data, provides, among other things, three 
sets of output. 

The potential implementations of the method, system, and computer program 
20 product of the invention is that it provides the following advantages, but are not limited 
thereto. First, the invention enhances existing home BG monitoring devices by 
producing and displaying: 1) estimated categories for HbAic, 2) estimated probability 
for SH in the subsequent six months, and 3) estimated short-term risk of hypoglycemia 
(i.e. for the next 24 hoxirs). The latter may include warnings, such as an alarm, that 
25 indicates imminent hypoglycemic episodes. These three components can also be 

integrated to provide continuous information about the glycemic control of individuals 
with diabetes, and to enhance the monitoring of their risk of hypoglycemia. 

As a second advantage, the invention enhances existing software or hardware 
that retrieves SMBG data. Such software or hardware is produced by virtually every 
30 manufacturer of home BG monitoring devices and is customarily used by patients and 
health care providers to interpret SMBG data. The methods and system of the 
invention can be directly incorporated into existing home blood glucose monitors, or 
used for the enhancement of software that retrieves SMBG data, by introducing a data 
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interpretation component capable of predicting both HbAu and periods of increased 
risk of hypoglycemia. 

Still yet another advantage, the invention evaluates the accuracy of home BG 
monitoring devices, both in the low and high BG ranges, and over the entire BG scale. 
5 Moreover, another advantage, the invention evaluates the effectiveness of 

various treatments for diabetes. 

Further still, as patients with diabetes face a life-long optimization problem of 
maintaining strict glycemic control without increasing their risk of hypoglycemia, the 
present invention alleviates this related problem by use of its simple and reliable 
10 methods, i.e., the invention is capable of evaluating both patients' glycemic control 
and their risk of hypoglycemia, and at the same time applying it in their everyday 
environments. 

Additionally, the invention provides the missing link by proposing three 
distinct, but compatible, algorithms for evaluating HbAic and the risk of hypoglycemia 

15 from SMBG data, to be used to predict the short-term and long-term risks of 
hypoglycemia, and the long-term risk of hyperglycemia. 

Finally, another advantage, the invention evaluates the effectiveness of new 
insulin or insulin delivery devices. Any manufacturer or researcher of insulin or 
insulin delivery devices can utilize the embodiments of the invention to test the relative 

20 success of proposed or tested insulin types or device delivery designs. 

The invention may be embodied in other specific forms without departing from 
the spirit or essential characteristics thereof. The foregoing embodiments are therefore 
to be considered in all respects illustrative rather than limiting of the invention 
described herein. Scope of the invention is thus indicated by the appended claims 

25 rather than by the foregoing description, and all changes which come within the 
meaning and range of equivalency of the claims are therefore intended to be embraced 
therein. 
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CLAIMS 
We claim: 

1 LA computerized method for evaluating the HbAic of a patient based on BG 

2 data collected over a predetermined duration, said method comprising; 

3 computing weighted deviation toward high blood glucose (WR) and estimated rate of 

4 change of blood glucose (Dr) based on said collected BG data; and 

5 estimating HbAic using a predetermined mathematical formula based on said 

6 computed WR and Dr. 

1 2. The method of claim 1, wherein: 

2 said computed WR is mathematically defined from a series of BG readings x/, x^, ... jc„ 

3 taken at time points tj^ t2, as: 

5 where: 

6 wr(BG;b) = 10,f(BG)^ if f(BG)>0 and 0 otherwise, 

7 b =1 , representing a weighting parameter, and 

8 

9 said computed DR is mathematically defined as: 

10 Dr = average of s^^-i - Sk, 

1 1 where: 

12 Sk^lO.S(k^~tif for k^OJ,..., t„-t J, 

13 S(tj)=f(xj), for j=l„.,n. 

1 3 . The method of claim 1 , wherein said the estimate of HbAic from said BG 

2 monitoring data is mathematically defined as : 

3 Estimated HbAic = 0.9008(WR) - 0.8207(DR) + 6.7489. 

1 4. The method of claim 1 , further comprising: 

2 defining predetermined categories for the estimate of HbAic , each of said HbAic 

3 estimate categories representing a range of values for estimated HbAic; and 

assigning said estimated HbAic to at least one of said HbAic estimate categories. 
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1 5. The method of claim 4, wherein said HbAic estimate categories are defined as 

2 follows: 

3 classified category 1, wherein said estimated HbAic is less than about 7.8; 

4 classified category 2, wherein said estimated HbAic is between about 7.8 and about 

5 8.5; 

6 classified category 3, wherein said estimated HbAic is between about 8.5 and about 

7 9.0; 

8 classified category 4, wherein said estimated HbAic is between about 9.0 and about 

9 9.6; 

10 classified category 5, wherein said estimated HbAic is between about 9.6 and about 

11 10.3; 

12 classified category 6, wherein said estimated HbAic is between about 10.3 and about 

13 11.0; and 

14 classified category 7, wherein said estimated HbAic is above about 1 1.0. 

1 6. The method of claim 5, further comprising: 

2 defining predicted confidence intervals for corresponding said HbAic estimate 

3 categories, wherein said predicted confidence intervals are defined as follows: 

4 said classified category 1 corresponds with a predicted HbAic less than about 

5 8.0; 

6 said classified category 2 corresponds with a predicted HbAic between about 

7 8.0 and about 8.5; 

8 said classified category 3 corresponds with a predicted HbAic between about 

9 8.5 and about 9.0; 

10 said classified category 4 corresponds with a predicted HbAic between about 

11 9.0 and about 9.5; 

12 said classified category 5 corresponds with a predicted HbAu between about 

13 9.5 and about 10.1; 

14 said classified category 6 corresponds with a predicted HbAic between about 

15 10.1 and about 11.0; and 
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1^ said classified category 7 corresponds with a predicted HbAic above about 

17 11.0. 

1 7. The method of claim 4, further comprising: 

2 defining predicted confidence intervals for corresponding said HbAic , each of said 

3 predicted confidence intervals representing a range of values for HbAic. 

1 8. The method of claim 7, wherein said predicted HbAic confidence intervals 

2 have about a 95% confidence level. 

1 9. A computerized method for evaluating the HbAic of a patient based on BG 

2 data collected over a predetermined duration, said method comprising: 

3 computing weighted deviation toward high blood glucose (WR) and estimated rate of 

4 change of blood glucose (Dr) based on said collected BG data; 

5 estimating HbAic using a predetermined mathematical formula based on said 

6 computed WR and Dr; and 

7 providing a predetermined confidence interval for classification of said estimated 

8 value of HbAic. 

1 10. The method of claim 9, wherein: 

2 said confidence interval is between about 85% to about 95%. 

1 1 1. A system for evaluating HbAic of a patient based on BG data collected over a 

2 predetermined duration, said system comprising: 

3 a database component operative to maintain a database identifying said BG data; 

4 a processor programmed to: 

5 compute weighted deviation toward high blood glucose (WR) and estimated 

6 rate of change of blood glucose (Dr) based on said collected BG data; and 

7 estimate HbAic using a predetermined mathematical formula based on said 

8 computed WR and Dr. 
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1 12. The system of claim 1 1 , wherein: 

2 said computed WR is mathematically defined from a series of BG readings x;, jc^, ... 

3 x,t taken at time points fy, t2, t„ as: 

5 where: 

6 wr(BG;b) = 10.f(BG/ if f(BG)>0 and 0 otherwise, 

7 b =1 , representing a weighting parameter, and 
8 

9 said computed DR is mathematically defined as: 

10 Dr = average of Sk+j - Sk, 

1 1 where: 

12 sk -^10.S(k+ tif for A:=ft 7, r„ - //, 

13 S(tj)=f(x]), for J=l n, 

1 13. The system of claim 1 1 , wherein said the estimate of HbAic from said BG 

2 monitoring data is mathematically defined as : 

3 Estimated HbAu = 0.9008(WR) - 0.8207(DR) + 6.7489. 

1 14. The system of claim 11, wherein said processor being further programmed to: 

2 define predetermined categories for the estimate of HbAic , each of said HbAic 

3 estimate categories representing a range of values for estimated HbAic; and 

4 assign said estimated HbAic to at least one of said HbAic estimate categories. 

1 15. The system of claim 14, wherein said HbAu estimate categories are defined as 

2 follows: 

3 classified category 1, wherein said estimated HbAic is less than about 7.8; 

4 classified category 2, wherein said estimated HbAic is between about 7.8 and about 

5 8.5; 

6 classified category 3, wherein said estimated HbAic is between about 8.5 and about 

7 9.0; 
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8 classified category 4, wherein said estimated HbAic is between about 9.0 and about 

9 9,6; 

10 classified category 5, wherein said estimated HbAic is between about 9.6 and about 

11 10,3; 

12 classified category 6, wherein said estimated HbAic is between about 10,3 and about 

13 11.0; and 

classified category 7, wherein said estimated HbAic is above about 1 1 .0. 

1 1 6. The system of claim 15, wherein said processor being further programmed to: 

2 define predicted confidence intervals for corresponding said HbAic estimate 

3 categories, wherein said predicted confidence intervals are defined as follows: 

4 said classified category 1 corresponds with a predicted HbAic less than about 

5 8.0; 

6 said classified category 2 corresponds with a predicted HbAic between about 

7 8.0 and about 8.5; 

8 said classified category 3 corresponds with a predicted HbAic between about 

9 8.5 and about 9.0; 

10 said classified category 4 corresponds with a predicted HbAu between about 

11 9.0 and about 9.5; 

12 said classified category 5 corresponds with a predicted HbAic between about 

13 9.5 and about 10.1; 

14 said classified category 6 corresponds with a predicted HbAic between about 

15 10.1 and about 11.0; and 

16 said classified category 7 corresponds with a predicted HbAic above about 

17 11.0. 

1 1 7. The system of claim 14, wherein said processor being further programmed to: 

2 define predicted confidence intervals for corresponding said HbAic » each of said 

3 predicted confidence intervals representing a range of values for HbAic. 
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1 18, The system of claim 1 7, wherein said predicted HbAic confidence intervals 

2 have about a 95% confidence level. 

1 1 9. A glycemic control system for evaluating HbAic of a patient, said system 

2 comprising: 

3 a BG acquisition mechanism, said acquisition mechanism configured to acquire BG 

4 data from the patient, 

5 a database component operative to maintain a database identifying said BG data; 

6 a processor programmed to: 

7 compute weighted deviation toward high blood glucose (WR) and estimated 

8 rate of change of blood glucose (Dr) based on said collected BG data; and 

9 estimate HbAic using a predetermined mathematical formula based on said 
10 computed WR and Dr. 

1 20. A computer program product comprising a computer useable medium having 

2 computer program logic for enabling at leeist one processor in a computer system to evaluate 

3 HbAic of a patient based on BG data, said computer program logic comprising: 

4 computing weighted deviation toward high blood glucose (WR) and estimated rate of 

5 change of blood glucose (Dr) based on said collected BG data; and 

6 estimating HbAic using a predetermined mathematical formula based on said 

7 computed WR and Dr. 

1 21. The computer program product of claim 20, wherein said computer program 

2 logic further comprises: 

3 providing a predetermined confidence interval for classification of said estimated 

4 value of HbAic, wherein said confidence interval is a single value or a range of values. 

1 22. A computerized method for evaluating the long term probability for severe 

2 hypoglycemia (SH) of a patient based on BG data collected over a predetermined duration, 

3 said method comprising: 
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4 computing weighted deviation toward low blood glucose (WL) and estimated rate of 

5 fall of blood glucose in the low BG range (DrDn) based on said collected BG data; and 

6 estimating the number of future SH episodes using a predetermined mathematical 

7 formula based on said computed WL and DrDn. 

1 23. The method of claim 22, wherein: 

2 said computed WL is mathematically defined from a series of BG readings xj, JC2, ... 

3 Xn taken at time points tj, t2, tn as: 

4 WL=^Twlixra) 

5 where: 

6 wl(BG;a) = 10J(BGf if f(BG)>0 and 0 otherwise, 

7 a =2, representing a weighting parameter, and 

8 said computed DR is mathematically defined as: 

9 DrDn = average of Sk^j - Sk, provided that Sk< Sk+i, 

10 where: Sk =10,S(k+ uf for A:=ft/, tn - tu 

11 S(tj)=f(Xj),ioxj=l...,n, 

1 24. The method of claim 22, wherein said estimated number of fiiture SH 

2 episodes (EstNSH) is mathematically defined as : 

3 EstNSH = 3.3613(WL) - 4.3427(DrDn) - 1.2716. 

1 25. The method of claim 22, fiirther comprising: 

2 defining predetermined EstNSH categories, each of said EstNSH categories 

3 representing a range of values for EstNSH; and 

4 assigning said EstNSH to at least one of said EstNSH categories. 

5 26. The method of claim 25, wherein said EstNSH categories are defined as 

6 follows: 

7 category 1, wherem said EstNSH category is less than about 0.775; 

8 category 2, wherein said EstNSH category is between about 0.775 and about 3.750 ; 
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9 category 3, wherein said EstNSH category is between about 3.750 and about 7.000; 

10 and 

1 1 category 4, wherein said EstNSH category is above about 7.0. 

1 27. The method of claim 26, further comprising: 

2 defining a probability of incurring a select number of SH episodes respectively for 

3 each of said assigned EstNSH categories; wherein said probability and said respective select 

4 number of SH are defined as: 

5 said classified category 1 corresponds with about a 90% probability of 

6 incurring about 0 SH episodes and about a 10% probability of incurring about 1 or 

7 more SH episodes over the predetermined duration; 

8 said classified category 2 corresponds with about a 50% probability of 

9 incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 
^ 0 episodes, and 25 % probability of incurring more than 2 SH episodes over the 

1 1 predetermined duration; 

12 said classified category 3 corresponds with about a 25% probability of 

13 incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 

14 episodes, and 50 % probability of incurring more than 2 SH episodes over the 

15 predetermined duration; and 

16 said classified category 4 corresponds with about a 20% probability of 

17 incurring about 0 to about 2 SH episodes and about a 80% probability of incurring 

18 more than 2 SH episodes over the predetermined duration. 

1 28. The method of claim 25, fiirther comprising: 

2 defining a probability of incurring a select number of SH episodes respectively for 

3 each of said assigned EstNSH categories; and 

4 providing at least one probability of incurring a select number of SH episodes 

5 according to said EstNSH category to which said EstNSH is assigned. 
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1 29. A computerized method for evaluating the long temi probability for severe 

2 hypoglycemia (SH) of a patient based on BG data collected over a predetermined duration, 

3 said method comprising: 

4 computing weighted deviation toward low blood glucose (WL) and estimated rate of 

5 fall of blood glucose in the low BG range (DrDn) based on said collected BG data; 

6 estimating the number of future SH episodes using a predetermined mathematical 

7 formula based on said computed WL and DrDn; and 

8 defining a probability of incurring a select number of SH episodes respective to said 

9 estimated SH episodes. 

1 30. A system for evaluating the long term probability for severe hypoglycemia 

2 (SH) of a patient based on BG data collected over a predetermined duration, said system 

3 comprising: 

4 a database component operative to maintain a database identifying said BG data; 

5 a processor programmed to: 

6 computing weighted deviation toward low blood glucose (WL) and estimated 

7 rate of fall of blood glucose in the low BG range (DrDn) based on said collected BG 

8 data; and 

9 estimating the number of future SH episodes using a predetermined 
10 mathematical formula based on said computed WL and DrDn. 

1 31. The system of claim 30, wherein: 

2 said computed WL is mathematically defined from a series of BG readings xj, X2, ... 

3 Xn taken at time points tj, t2, tn as: 

4 WL=^Y.wliixra) 

5 where: 

6 wl(BG;a) = 10.f(BGf if f(BG)>0 and 0 otherwise, 

7 a =2, representing a weighting parameter, and 
8 
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9 said computed DR is mathematically defined as: 

10 DrDn = average of Sk^i - Sk, provided that Sk< 5jt+/, 

1 1 where: 

12 Sk =10,S(k+ tif for k^OJ, tn - tu 

13 5r/y;==/rxy),fory=7,...,/2. 

1 32. The system of claim 30, wherein said estimated nvimber of future SH episodes 

2 (EstNSH) is mathematically defined as : 

3 EstNSH = 3.3613(WL) - 4.3427(DrDn) - 1.2716. 

1 33. The system of claim 30, wherein said processor being further programmed to: 

2 define predetermined EstNSH categories, each of said EstNSH categories 

3 representing a range of values for EstNSH; and 

4 assign said EstNSH to at least one of said EstNSH categories. 

1 34. The system of claim 33, wherein said EstNSH categories are defined as 

2 follows: 

3 category 1, wherein said EstNSH category is less than about 0.775; 

4 category 2, wherein said EstNSH category is between about 0.775 and about 3.750 ; 

5 category 3, wherein said EstNSH category is between about 3.750 and about 7.000; 

6 and 

7 category 4, wherein said EstNSH category is above about 7.0. 

1 35. The method of claim 34, wherein said processor being further programmed to: 

2 define a probability of incurring a select number of SH episodes respectively for each 

3 of said assigned EstNSH categories, wherein said probability and said respective select 

4 number of SH are defined as: 

5 said classified category 1 corresponds with about a 90% probability of 

6 incurring about 0 SH episodes and about a 10% probability of incurring about 1 or 

7 more SH episodes over the predetermined duration; 
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8 said classified category 2 corresponds with about a 50% probability of 

9 incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 

10 episodes, and 25 % probability of incurring more than 2 SH episodes over the 

1 1 predetermined duration; 

12 said classified category 3 corresponds with about a 25% probability of 

13 incurring about 0 SH episodes, 25 % probability of incurring about 1 to about 2 SH 

14 episodes, and 50 % probability of incurring more than 2 SH episodes over the 

15 predetermined duration; and 

16 said classified category 4 corresponds with about a 20% probability of 

17 incurring about 0 to about 2 SH episodes and about a 80% probability of incurring 

1 8 more than 2 SH episodes over the predetermined duration. 

1 36. The system of cleiim 33, wherein said processor being further programmed to: 

2 define a probability of incurring a select number of SH episodes respectively for each 

3 of said assigned EstNSH categories; and 

4 provide at least one probability of incurring a select number of SH episodes according 

5 to said EstNSH category to which said EstNSH is assigned. 

1 37. A glycemic control system for evaluating the long term probability for severe 

2 hypoglycemia (SH) of a patient, said system comprising: 

3 a BG acquisition mechanism, said acquisition mechanism configured to acquire BG 

4 data from the patient, 

5 a database component operative to maintain a database identifying said BG data; 

6 a processor programmed to: 

7 computing weighted deviation toward low blood glucose (WL) and estimated 

8 rate of fall of blood glucose in the low BG range (DrDn) based on said collected BG 

9 data; and 

10 estimating the number of future SH episodes using a predetermined 

1 1 mathematical formula based on said cpmputed WL and DrDn. 
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1 38. A computer program product comprising a computer useable medium having 

2 computer program logic for enabling at least one processor in a computer system to evaluate 

3 long term probability for severe hypoglycemia (SH) of a patient based on BG data, said 

4 computer program logic comprising: 

5 computing weighted deviation toward low blood glucose (WL) and estimated rate of 

6 fall of blood glucose in the low BG range (DrDn) based on said collected BG data; and 

7 estimating the number of future SH episodes using a predetermined mathematical 

8 formula based on said computed WL and DrDn. 

1 39. The computer program product of claim 38, wherein said computer program 

2 logic further comprises: 

3 defining a probability of incurring a select mmiber of SH episodes respective to said 

4 estimated SH episodes. 

1 40. A computerized method for evaluating the short term risk for severe 

2 hypoglycemia (SH) of a patient based on BG data collected over a predetermined duration, 

3 said method comprising: 

4 computing weighted deviation toward low blood glucose (WL); 

5 determining Max(wl) by calculating maximum value of wl(BG;2); and 

6 determining risk value by taking the geometric mean of WL and Max(wl) over said 

7 predetermined duration, said risk value is mathematically defined as: 

8 risk value = -^WL • Max(wl) . 

1 41. The method of claim 40, wherein: 

2 said computed WL is mathematically defined from a series of BG readings xj, X2, ... jc„ 

3 taken over the predetermined duration as: 

4 wL^^Ys^iixra) 

5 where: 

6 wl(BG;a)= 10/(BG/ if /(«G;>0 and 0 otherwise, 

7 a =2, representing a weighting parameter. 
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1 42. The method of claim 40, further comprising: 

2 providing a predetermined threshold risk value; and 

3 comparing said determined risk value to said threshold risk value. 

1 43 . The method of claim 42, wherein: 

2 if said determined risk value is greater than said threshold value then short term risk 

3 of incurring a hypoglycemic episode is high; and 

4 if said determined risk value is less than said threshold value then short term risk of 

5 incurring a hypoglycemic episode is low. 

1 44. The method of clmm 43, wherein said short term is approximately a 24 hour 

2 period. 

1 45. The method of claim 43, wherein said short term ranges from about 12 to 

2 about 72 hour period. 

1 46. The method of claim 43, wherein said threshold value is approximately 17. 

1 47. The method of claim 43, wherein said threshold value is between about 12 to 

2 25. 

1 48. A system for evaluating the short term risk for severe hypoglycemia (SH) of a 

2 patient based on BG data collected over a predetermined duration, said system comprising: 

3 a database component operative to maintain a database identifying said BG data; 

4 a processor programmed to : 

5 compute weighted deviation toward low blood glucose (WL); 

6 determine Max(wl) by calculating maximirai value of wl(BG;2)\ and 

7 determine risk value by taking the geometric mean of WL and Max(wl) over 

8 said predetermined duration, said risk value is mathematically defined as: 

9 risk value = 
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1 49. The system of claim 48, wherein: 

2 said computed WL is mathematically defined from a series of BG readings jc/, xp, ... x„ 

3 taken over the predetermined duration as: 

5 where: 

6 wl(BG;a) = 10.f(BGf if f(BG)>0 and 0 otherwise, 

7 a =2, representing a weighting parameter. 

1 50. The system of claim 48, wherein said processor being further programmed to: 

2 provide a predetermined threshold risk value; and 

3 compare said determined risk value to said threshold risk value. 

1 51. The system of claim 50, wherein: 

2 if said determined risk value is greater than said threshold value then short term risk 

3 of incurring a hypoglycemic episode is high; and 

4 if said determined risk value is less than said threshold value then short term risk of 

5 incurring a hypoglycemic episode is low. 

1 52. The system of claim 51, wherein said short term is approximately a 24 hour 

2 period. 

1 53. The method of claim 51, wherein said short term ranges from about 12 to 

2 about 72 hour period. 

1 54. The system of claim 5 1 , wherein said threshold value is approximately 1 7. 

1 55 . The method of claim 5 1 , wherein said threshold value is between about 1 2 to 

2 25. 
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1 56. A glycemic control system for evaluating the short term risk for severe 

2 hypoglycemia (SH) of a patient, said system comprising: 

3 a BG acquisition mechanism, said acquisition mechanism configured to acquire BG 

4 data from the patient, 

5 a database component operative to maintain a database identifying said BG data; 

6 a processor programmed to: 

7 compute weighted deviation toward low blood glucose (WL); 

8 determine Max(wl) by calculating maximxmi value of wl(BG;2); and 

9 determine risk value by taking the geometric mean of WL and Max(wl) over 

10 said predetermined duration, said risk value is mathematically defined as: 

1 1 risk value = -yJJVL • Max(wl) . 

1 57, A computer program product comprising a computer useable medium having 

2 computer program logic for enabling at least one processor in a computer system to evaluate 

3 the short term risk for severe hypoglycemia (SH) of a patient based on BG data collected over 

4 a predetermined duration, said computer program logic comprising: 

5 computing weighted deviation toward low blood glucose (WL); 

6 determining Max(wl) by calculating maximum value of wl(BG;2); and 

7 determining risk value by taking the geometric mean of WL and Max(wl) over said 

8 predetermined duration, said risk value is mathematically defined as: 

9 risk value = -yJWL • Max(wl) . 

1 58. The computer program product of claim 57, wherein said computer program 

2 logic further comprises: 

3 providing a predetermined threshold risk value; and 

4 comparing said determined risk value to said threshold risk value. 
5 
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